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ABSTRACT 

We investigate the decay mechanisms in the D"^ — > wk^ and -D^ p°7r+ 
transitions. The naive factorization ansatz predicts vanishing amphtude for 
the ujTi^ decay, while the p^ir^ decay amphtude does have 

an annihilation contribution also in this limit. Both decays can proceed 
through intermediate states of hidden strangeness, e.g. K,K*, which we 
estimate in this paper. These contributions can explain the experimental 
value for the un^ decay rate, which no longer can be viewed as a 

clean signature of the annihilation decay of Df. The combination of the 
7r(1300) pole dominated annihilation contribution and the internal K,K* 
exchange can saturate present experimental upper bound on p^n'^ 
decay rate, which is therefore expected to be within the experimental reach. 
Finally, the proposed mechanism of hidden strangeness FSI constitutes only 
a small correction to the Cabibbo allowed decay rates Dg KK*, (pTT, which 
are well described already in the factorization approximation. 



1 Introduction 



It has been suggested [1] , that the observation of the —>■ wtt"*" decay 

^a;7r+) = (2.8± 1.1) X 10-^ (1) 

can be seen as a clean signature of the annihilation decay of Df. The sizes of 
these contributions are of great phenomenological interest, but are very hard 
to obtain from theoretical considerations. In particular, the factorization ap- 
proximation gives vanishing prediction for BR{D^ as an immediate 
consequence of the fact that due to the isospin and G parity = of cjtt"'" 
state there are no annihilation contributions with one intermediate resonant 
state (cf. also Fig. 1). The interpretation of branching ratio (1) as an (al- 
most) purely annihilation decay, would then also give important information 
about the, otherwise relatively poorly known, sizes of annihilation diagrams 
in other related hadronic decays. In this paper, however, we will argue that 
the experimental value for the U7r~^ transition can be accommodated 
(already) by considering only color suppressed spectator decay with subse- 
quent final state interactions (FSI). This leaves little room for unambiguous 
study of annihilation effects from the cjvr"'" decay mode. 

Another issue that we address is the size of the p^Ti^ transition. 

The experimental situation regarding this decay is somewhat unclear. The 
difficulties arise from the fact that the — > p^TT"*" decay is observed as a 
resonance in the three body decay — > 7r~7r"'"7r"'", which is dominated by the 
decays through isoscalar resonances, /o(980)7r"'", /o(1370)7r+ [2, 3]. The decay 
p^TT^ therefore constitutes just a small correction to these dominant 
processes. The two available experimental analyses give 

BR{Dt ^ p°7r+) < 7 X 10"^ [2] , (2) 

BR{D+ p\+) = (5.9 ± 4.6) x 10"^ [3] , (3) 

where the first is a 90% CL limit, while in the second case we have added all 
the errors in quadrature. Clearly, the two measurements are in agreement 
with each other. The error on the number from [3] is still very large and can 
at best be viewed as an indication toward possible value for BR{Df p^7i~^). 
In particular, it is still compatible with zero. As we will show in the following, 
however, the expectations for BR{Df — p°7r+) one gets from spectator 
decays with FSI and/or annihilation diagrams are just in the same ballpark 
as the above experimental indications. 

An important difference between ujti^ and p^n^ transitions 

are the the quantum numbers of the final state. The (utt^) is a 
1"'"(0^) state, while the (p°7r+) can be either in the I'^ — 1~ or 2~ state. 
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Figure 1: Annihilation diagram of D+ p°Tr+ decay with "it" = tt, 7r(1300), 7r(1800), 
with dot representing the weak vertex. 

again with JP = Q-. A scan through PDG book [4] reveals, that there are 
no resonant states with I'-^{J^) = l^(O^) assignment, while there are three 
known states with I^{J^) = 1^(0^), the two strongly decaying resonances 
7r(1300) and 7r(1800), and the pion. 

The presence of resonances near the Dg meson mass with the quantum 
numbers of (p°7r^) state indicates an enhancement of annihilation contri- 
butions, which in this case can proceed through an intermediate resonating 
state (see Fig. 1). To gain an insight on the numerical importance of the 
annihilation processes we estimate the contribution coming from the 7r(1300) 
intermediate state. To do so, we need to estimate the decay constant of 
7r(1300). Taking the PDG [4] upper bound for r 7r(1300)z/^ one arrives 
at /7r(i3oo) < 4 MeV. In the factorization approximation for the weak vertex 
(for more details see section 2) we then get 

BR{Dt ^ p°7r+),(i3oo) < 7 x 10"^ , (4) 

where we have used fo^ = 230 MeV, together with the conservative as- 
sumptions of i?-R(7r(1300) pn) ~ 100% and r(7r(1300)) equal to its upper 
experimental bound of 600 MeV. Most probably this slightly overestimates 
the contribution of 7r(1300) intermediate resonant state to the p°7r+ 
decay width, as also the presence (but not the size) of other decay channels 
of 7r(1300) has been seen experimentally (e.g. 7r(1300) {7ni)s-waye'^ [4]). 
Note in particular that the above assumptions about the 7r(1300) — >■ pvr de- 
cay width correspond to the case where 5'p7r7r(i3oo) — gpTrn- We therefore do 
not expect the 7r(1300) contribution to BR{Df p°7r+) to lie significantly 
below the upper limit (4). 

The importance of the upper limit (4) is that the contribution from 
7r(1300) can even saturate the 90 % CL experimental upper bound (2). Of 
course the actual size of 7r(1300) contribution is not known and lies some- 
where below the upper bound (4). Also, interference with other annihilation 
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contributions from intermediate tt and 7r(1800) states can somewhat change 
the above estimate (using PCAC, the contribution from tt was found in [5] 
to be negUgible, while the contribution of 7r(1800) is difficult to estimate due 
to the lack of experimental data). Furthermore, as we will show in the next 
section, the contributions of final state interactions fall in exactly the same 
range (cf. Eq. (15) below). The lesson to be learned from this simple exercise 
is that, unless there are large cancellations, the value of BR{Df — > p^TT"*") is 
expected to be near to its present experimental upper bound and should be 
measured by the FOCUS collaboration in the near future [6]. 

On the other hand note that, as explained above, no such resonance 
enhancements of annihilation contributions are possible for the (cutt"*") fi- 
nal state. How can then one explain a relatively large experimental value 
for BR{Ds ^7r+), Eq. (1)? The answer lies in the fact that there are 
multi-body intermediate states that do have correct values of and J^, for 
instance the two-body {K^*^K^*^) states (see Fig. 2). Moreover, such states 
of hidden strangeness can be obtained from a (color suppressed) spectator 
decay of the Dg meson and are therefore expected to be sizable. We will 
estimate the sizes of these contributions in the next section. 



2 Hidden strangeness FSI 

In estimating the contributions from hidden strangeness intermediate states 
that can arise from spectator quark diagrams, we will resort to the following 
simplifications 

• Only two body intermediate states with s, s quantum numbers will be 
taken into account. Moreover, only the contributions of lowest lying 

pseudoscalar and vector states (neglecting their decay widths) will be 
considered. Note that the re-scattering through intermediate K, K* 
states is possible for both p^TT"*" as well as uj-n^ final state (cf. Fig. 2), 
while the re-scattering with intermediate r/ or rj' is possible only in the 
case of uj-K^ final state due to isospin and G parity conservation (cf. 
Fig. 3). 

• For the weak transition L>+ ^ (irWxW)+ intheL'+ ^ (XWxW)+ ^ 
p°7r+ and Df {K^*^ K^*'>)~^ cj7r+ decay chains as well as for 
the weak transition ri{ri')p^ in the ri{ri')p^ cjtt"'" 
decay chain we will use the factorization approximation. The weak 
Lagrangian is therefore 

>Cweak = — -i=VcsV*a{ai{ud)H{sc)H + a2{sd) h{uc) h) , (5) 
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Figure 3: The intermediate 77, r]' contributions in the Z)+ — > a;7r+ decay. 
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decay 




BRth 


Dt K*+K^* 


(5.8 ±2.5)% 


6.4% 


Dt K+K"* 


(3.3 ±0.9)% 


3.4% 


Dt K*+K^ 


(4.3 ±1.4)% 


2.7% 


Dt VP^ 


(10.8 ±3.1)% 


5.6% 


Dt ^ V'P^ 


(10.1 ±2.8)% 


2.2% 



Table 1: The branching ratios for the Z)+ two body decays: in the second column 
the experimental results are listed, while in the third the predictions are given when the 
factorization approximation is used. 



with {ud)H, ■ ■ ■ the hadronized V-A weak currents, Vij the CKM matrix 
elements and ai^2 the effective (phenomenological) Wilson coefficients 
taken to be ai — 1.26 and 02 = —0.52 [7-9]. 

• Finally, the strong interactions are taken into account through the fol- 
lowing effective Lagrangian [10-12]: 



Strong = ^Tr(p'^[n, d,U]) - A^e'^'^^^Trid.p^d^ppU) , (6) 

where 11 and are 3x3 matrices containing pseudoscalar and vector 
meson operators respectively and / is a pseudoscalar decay constant 
(11). We use numerical values Cyvu = 0.33, and gp-^^j- = 5.9 [10-12]. 

We have checked that the use of factorization (at tree level, with values 
of form factors as given below in Eqs. (7)-(ll) and in Table 2) for the 
Dt K*+K*^, Dt K+K*^ and ^ K^K*+ decays gives reasonable 
estimates of the measured rates (note that we do not need Dt — > K'^K'^ in 
further considerations), with the results listed in Table 1. In these results 
the annihilation contributions have been neglected since they are an order of 
magnitude smaller. We also neglect the triangle graphs of the sort shown on 
Figs. 2, 3 as we expect them (based on the numerical results for u;7r~^, p^n^) 
to be suppressed relative to the tree level contributions. 

The situation in the case of 77, 77' intermediate states is not so favorable. 
To treat the r], rj' mixing we use the approach of ref. [13] with the value of 
the mixing angle transforming between 77,77' and 77^ ~ {uu + dd)/\/2,r]s ~ 
ss states taken to be ^ = 40°. The factorization approach then gives a 
reasonable description of Dt — > ^"''77 decay, while it does not reproduce 
satisfactorily the experimental result for Dt p'^rj' (cf. Table 1). This 
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is a known problem as Df — > rate is very difficult to reproduce by 
any of present approaches [8,9, 14]. This inevitably introduces some further 
uncertainty into our approach, yet the resulting uncertainty is not expected 
to affect significantly our main conclusions. 

For the weak current matrix elements between Dg and vector or pseu- 
doscalar final states we use a common decomposition 



mp-\-mv q 
+i{M + my) \e^-^-^q^] A,{q') - i \p^^ - ^^-ZJ^q^] A^{q^) , 



g2 



q^ 



(7) 

(8) 

where is a vector meson characterized by the polarization vector and 
mass my, while pseudoscalar mesons P, Dg have masses m, M. We use fur- 
ther abbreviations F'" = 7^'(1 — 7^), qf^ = — k^^ and P^ = p^ + . 

For the q^ dependence of the form factors we use results of [15], based on 
a quark model calculation combined with a fit to lattice and experimental 
data. Ref. [15] provides a simple fit to their numerical results with the form 
factors V{q^) and Ao^q^) described by double pole q^ dependence 

' (1 - gVM2)(l - ' ^ > 

while single pole parameterization 

can be used for Ai^2{q'^), as the contributing resonance lie farther away from 
the physical region (note that this parameterization applies also to Fq form 
factor, which however does not contribute in the processes we discuss in this 
paper). The values of /(O) and a are fisted in Table 2 and are taken from 
[15]. We use M = 1.97 GeV in the expression for Aq, while M = 2.11 GeV is 
used for other form factors [15]. Incidentally, the parameterizations of form 
factors (9), (10) make all the loop diagrams in Figs. 2 and 3 finite, so that no 
regularization procedure is needed (we have also checked that the numerical 
results do not change significantly, if one uses a single pole parameterization 
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fonii factor 




V 


-4o 


-4i 


-42 




/(O) 


0.72 


1.04 


0.67 


0.57 


0.42 


0.78 


a 


0.2 


0.24 


0.2 


0.29 


0.58 


0.23 



Table 2: Form factors at = [15]. The results in the first five columns are for 
Ds K, K*lvi transitions. The last column stands for the form factor appearing in 
Ds rjslvu (the ss component oiri,ri') transition. 



D+ un^ 


Ai 


^2 


A, 


Aa 


A, 


A 




Bfj' 


•AiD 


-0.7 


0.7 


-1.1 


-1.4 


11.3 


12.5 


1.3 


3.6 


AiA 


-0.7 


0.7 


3.3 


1.5 


-4.0 


-19.7 


-7.2 


-3.7 



Table 3: The dispersive AiD and absorptive AiA parts of the amplitudes (in units of 
10~^ GeV) for the Df — > wtt^ decay corresponding to the diagrams on Fig. 2 (Ai) and 
Fig. 3 {Brj.T]')- The amplitudes for the Df /3°7r+ decay (neglecting the mass difference 
between rup and moj) are obtained by inverting the sign of Aid, AiA for even i, while 
^ri,ri' = 0. 



instead of (9) and then uses a cut-off regularization with some scale above 
but close enough to meson mass) . 

For the decay constants, defined through 

mr^,q\P{p)) = ifpp" , {0\qrq\V{p)) = gve^ ■ (11) 

we use the following values: fo — 0.207 GeV and fos = 1-13/d as obtained 
on the lattice [16] and for the rest fx = 0.16 GeV, \gK*\ = 0.19 GeV^ 
\gp\ — 0.17 GeV^ and \g^\ — 0.15 GeV^ coming from the experimental mea- 
surements [4]. 

The amplitudes for the ujt:^ and Df —>■ p^TT"*" decays can be written 

as: 

k2{J2A^^ + B), (12) 

i 

k,Y.4'\ (13) 

i 

with e the helicity zero polarization vector of the uj or p vector mesons, 
while k2 is the pion momentum. The reduced amplitudes A^'^'^'^^^ and B 
correspond to the diagrams on Figs. 2 and 3 respectively (note that in our 
phase convention for Clebsch-Gordan coefficients Ai^^ = (— l)*"'"^^-'^^ i.e. the 
even diagrams for the two channels on Fig. 2 differ by a sign) . The explicit 



A{Dtip)^u;{e,h)7r+{k2)) 



Gi 



A{Dt{p)^p\eMW{k2))^^e- 
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expressions can be found in Appendix. The numerical values for Af and 
B are given in Tabic 3. 

Combining the above results we arrive at the prediction 

BR{Dt u;7r+) = 3.0 x 10"^ (14) 

Note that in this calculation we have used the factorization approximation for 
the diagram of Fig. 3, which as stated above, does not work well for —>■ 
p^r],r]' transition.^ If one instead uses experimental input to rescale the 
corresponding amplitudes one ends up with BR(Df — > a;7r+) = 4.4 x 10~^. 
We see, that it is actually very easy to explain the experimental result (1) 
entirely in terms of spectator quark transition together with FSl. The claim, 
that D'^ — > ujt:~^ can be used as a probe of annihilation contributions in 
hadronic physics is therefore not justified. 

As what concerns the Df — > p°7r"^ transition, the FSI contributions alone 
give 

BR{D+ p°7r+)Fsi = 0.7 x lO'^ , (15) 

which is almost exactly the same as our estimate of the upper bound on the 
annihilation contribution (4) and actually coincides with the present 90% 
CL upper bound. If there is no destructive interference between these two 
contributions and the contributions of FSI through higher resonances that we 
did not take into account, this decay mode should be established in the near 
future. This expectation is supported also by other theoretical approaches 
which give the rate for p^TT"*" to be equal [9] or even larger than the 

rate for —>■ lutt'^ decay [14, 17]. 

Possible cancellation that may occur, however, make the theoretical pre- 
dictions rather uncertain. Adding the FSI contribution and the maximal 
annihilation contributions (4) with alternating signs gives a fairly large in- 
terval 

BR{D+ p°7r+) = (0.05 - 3.5) x 10"^ . (16) 

We note that the experimental uncertainties translating in the input param- 
eters can change the values for BR{D'^ — > p^tt~^) and BR{Df — > u!7i~^) by 
about 20%. 

Finally, we mention that the kind of FSI contributions we were considering 
in this paper will not be the leading contribution in the (pn^ transition, 

namely, this decay can proceed through spectator quark transition directly. 
Use of factorization approximation for the weak vertex leads to a prediction 
BR{Df — > (/)7r+) = 4.0%, which is already in excellent agreement with the 

^Adding hidden strangeness FSI in the spirit of this paper to the Df p'^rj' channel 
does not mend the discrepancy, as it gives an order of magnitude smaller contribution. 
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experimental result of 3.6 ± 0.9%. Inclusion of FSI reduces the theoretical 
prediction from 4% to ~ 3.6% and does not spoil the agreement with the 
experiment (it actually even improves it). The size of the shift also indicates 
that FSI of the type described in the present paper are in the case of 
(j)-K^ transition a second order effect. Note as well, that the size of the FSI 
correction is in agreement with the predictions for BR{D~^ — > p°7r+) and 
BR{D~^ — > a;7r+), which are an order of magnitude smaller than BR{Df — > 

4)71+ ) . 

In conclusion, we found that the hidden strangeness final state interac- 
tions are very important in understanding the D+ — > ujn+ and D+ — > p^7i+ 
decay mechanism. The ujt:+ amplitude can be explained fully by this 

mechanism. As for the D+ p^Ti+ amplitude the predictions we obtain lie 
in a fairly large range due to possible cancellation between FSI and single 
pole contributions. There is a hope that D+ p^7r+ will be measured soon, 
hopefully shedding more light on our understanding of the D+ — > p^7i+ de- 
cay mechanism. Finally, we note that the hidden strangeness FSI discussed 
in this paper represents a second order effect, the inclusion of which does 
not spoil the good agreement of factorization approximation obtained for 
Ds (pn, KK*. 
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A Decay amplitudes 

In this appendix we list expressions for the reduced amplitudes Ai and B 
defined through Eqs. (12)-(13) and Figs. 2, 3. We present explicitly the 
Af'^ reduced amplitudes only for the D+ uj7t+ decay. The corresponding 
expressions A^'^^ for the Df — > p^7r+ decay are then obtained by inverting 
the signs of Af^^ for even i: 

r k A^-^ - I an ^ r "^(g""' " ^2^^ MP 
-J (2^"^^^m. + M1- pl/Mf'^-'^'^' 

p\-ml q^-m^ ^ /^V2 

(17) 
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p\ -ml q^-m^ ^ /^^ 

(18) 

J (27r)4 1 — P2/Mp pi — mi pf — 



-I 



(19) 



X 



PI / / p{-w? pi- ml 

-iig""' - g"g"7m^) i4Cyyn e^/g,/^ p^q^ iACyyu e^arS q^e^kf 
(f" - ml U V^fK 

(20) 

(^) /• A i^+(0) , ^ . -i(g^'--p^p^/m^) ^^ 

^ i igp7T7r{k2 - q)v igpTTTrjq - Pi) • £ 

pI — m^ q^ — m^ ^ 2 ' 



^ -J (27r)^"^ r P? 1 - Pl/Ml + 1 - P?/M| 

P+P2^ ^Pl {-tfKPl)x 2 2 X 

PI / / P2 - ^* 

i « ^fi'p7r7r(g - k2)u igpwAi - Pi) • £ 



(22) 

In expressions above we used m, m*, M for the K, K*, Ds mass respectively. 
The pole masses used in the form factors are Mp = 2.11 GeV and Ma = 1.97 
GeV. The definitions of momenta ki^2, Pi, 2 and q are the same as the ones 
used on Figs. 2, 3. 
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The B part of the Df un^ amphtude (Fig. 3) is 

, „ /"A ^-t..,<"(°) , ^ , -'{3'"'' -pVi) ' .. 



a; 

4^ (T T (5 

/»,(') 



(23) 



while ;B = for —>■ p^i:^ . In (23) we used the abbreviations 

Kr^ = cos 0/^2 X^/ = sin 0/^2 (24) 

and 

F+,,(0) = sin0F,^,+ (O) F+,,,(0) = cos0F,^,^(O). 
Finally, the decay rate for the D'^ — > a;(p°)7r''" is given by 



(25) 



with M the Dg mass, and mi^2 the vector and the pseudoscalar meson masses 
respectively, while 

k = ^[(M^ - (mi + m2)2)(M2 - (mi - m^f)]"' (26) 
is the size of the three-momentum of the final particles in the Dg rest frame. 
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